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Abstract 

The signals in high mass dilepton spectroscopy with nucleus-nucleus colli- 
sions at collider energies are superimposed on a generally large combinatorial 
background. Because this background contains a significant correlated like- 
sign component originating from B meson decays, the "like-sign" method to 
determine the background is inappropriate. We discuss strategies to deal with 
the correlations in the background. By taking advantage of the B meson os- 
cillation mechanism and of the particular features of B meson decays a new 
method to measure the bb production cross-section is proposed. 
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I. INTRODUCTION 



Heavy quarkonia states provide, via their leptonic decays, an essential tool to probe 
the earliest stage of heavy ion collisions at ultra-relativistic energies and are considered a 
crucial signature for demonstrating the existence of the quark gluon plasma (QGP) and for 
diagnosing its properties f^-^. The anomalous suppression of the J /if) meson, observed in 
Pb + Pb reactions at CERN/SPS, is indeed considered as an indication for the creation of 
a deconfined medium However, the production mechanism for both open and hidden 
charm at SPS is a subject of intensive current debate. Thus the standard pQCD-based 
framework [|l],^ for charmonium production and interaction with the deconfined medium is 
nowadays confronted with new approaches such as the statistical hadronization model |^ 
and the kinetic recombination model [0. These statistical models differ from the standard 
approach in their predictions for higher beam energies: they imply a significant J/ip en- 
hancement rather than suppression. In both approaches rescattering of heavy quarks in a 
deconfined medium is essential. Furthermore, a unique signature of the new mechanisms is 
the pecuhar centrality dependence 0, a scaling of the quarkonia yield with the square of 
the open charm (or bottom) yield and thermal ratios for relative yields of quarkonia 

states 0. It then becomes mandatory to measure, in addition to the quarkonia yields, also 
the yields for open charm and open beauty. The study of high mass dileptons is therefore 
one of the major physics goals of RHIC and the LHC heavy ion program. This new regime 
of beam energies brings however the challenge of extracting, for the first time, the quarkonia 
signals in the presence of a significant and highly complex combinatorial background which 
arises mostly from the semi-leptonic decay of open charm and open bottom (up to 200 cc 
and 6 bb pairs are expected to be produced per central Pb + Pb collision at LHC energy). 
Consequently, high mass dilepton spectra exhibit peculiar features and the well known tech- 
niques for background subtraction which work successfully for low mass resonances at low 
beam energies, cannot be applied in a straightforward way. 

Based on simulations performed for Pb + Pb collisions at LHC energy we discuss the 
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problems related to the subtraction of the combinatorial background to the dilepton spectra 
at high invariant mass. We demonstrate that the peculiar characteristics of the combinatorial 
background at high invariant mass provides a means to measure, albeit indirectly, the number 
of bb pairs produced in the reaction. 

II. SUBTRACTING THE COMBINATORIAL BACKGROUND 

In the multiparticle environment characteristic of heavy ion reactions the possibility to 
observe a resonance relies on the fact that the reconstructed invariant mass of its decay 
products appears as a narrow peak signal superimposed on a broad background. Depending 
on the underlying physics and on the event multiplicity, the background originates from 
uncorrelated particles and/or from correlated particles i.e. of common origin. In principle 
the signal can be extracted by fitting the signal+background distribution with appropriate 
functions chosen in order to provide a good description of the overall spectrum. This tech- 
nique, however, does not work if signal and background have a similar shape. This could 
result either i) from a bump in the background due to physical reasons, detector acceptance, 
or analysis cuts, or ii) from a broad signal due to a large intrinsic resonance width and/or 
a smearing of this width by the apparatus resolution. In these cases it becomes difficult or 
even impossible to disentangle the signal and background components by fitting their sum. 
Obviously, the situation gets even worse when the resonant signal is small and/or when the 
shape of the background is unknown. These difficulties become even more obvious when 
extracting a continuum like the Drell-Yan, the thermal radiation, or the open charm/bottom 
dilepton signal. The only way to overcome this problem is to estimate independently the 
background distribution and to subtract it from the signal+background spectrum. Such a 
technique is nowadays intensively used in heavy ion collisions at high beam energies. 

Considering an unlike-sign two particle invariant mass spectrum, there are mainly two 
ways to determine and subtract the combinatorial background: 

• Like-sign pair technique: the uncorrelated background in the unlike-sign pair sample 
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is estimated by the number of like-sign pairs within each event. The signal S of the 
number of unlike-sign pairs N"*"" is then given by S = N+- - 2RVN++N~ where the 
factor R accounts for a possible asymmetry in the production of positively and nega- 
tively charged particle and/or an asymmetry due to a detector trigger/acceptance bias 
relative to the particle charge. This technique has the disadvantage that the statis- 
tics in the background spectrum is limited to the number of available events. On the 
other hand, since the number of unlike-sign pairs and like-sign pairs are calculated 
within the same event, the normalization of the determined background to the sig- 
nal+background spectrum is straightforward, provided that one knows the R factor 
with good accuracy. 

• Event-mixing technique: the uncorrelated background in the unlike-sign pair sample 
is given by an unlike-sign pair spectrum for which the two particles of a pair are 
taken from two different events. This offers a better statistical precision than the 
previous method since one can mix each event with several (or even many) other 
events. The mixed events have to be identical in centrality. In presence of flow 
they have to be rotated into a common reaction plane. After normalization of the 
signal+background spectrum and the event-mixing spectrum to the respective number 
of events an additional normalization factor of 2 has to be applied to the event-mixing 
spectrum. Note that with real data, the normalization factor is not straightforward 
because, due to fluctuations and detector resolution, the centrality of the mixed events 
is never strictly identical. Therefore the normalization of the estimated background 
to the signal+background is usually done from the integration of the two spectra in a 
region where the correlated signal is assumed to be negligible. Note also that, in order 
to be applied to experimental data, the event-mixing technique requires that the two 
track resolution of the detector be taken into account. 

Each method has its own advantages and drawbacks depending on the particle environment 
and on the nature of the signal to be extracted [^HHI- In particular, the event-mixing 

4 



technique can lead to incorrect results in the case of large signal-to-background. This can 
be intuitively understood by considering the extreme case of a sample of events where each 
event consists of two correlated particles only, such as a cc pair. In this case there would 
be no real background to the signal but the event-mixing technique would produce a fake 
background. On the other hand it is obvious that the like-sign method, which relies on 
the fact that a like-sign pair is always uncorrelated, will fail if the events contain like-sign 
correlated particle pairs. 

In the following we shall show that, because the decay products of bottom decay contain 
like-sign correlated lepton pairs, the applicability of the like-sign technique for background 
subtraction in dilepton physics at the collider energies, such as currently available at RHIC 
and to become available at the LHC, is questionable. By making use of the differences 
between the background estimated with the like-sign and the event-mixing techniques, we 
propose a new method to measure indirectly the number of bb pairs from the dilepton 
spectra. As illustrative examples, we consider the two dilepton channels (e"'"e~, fi~^fJ'~) in 
the acceptance of the ALICE detector |T2[ for central Pb + Pb collisions. Technical details 
about the simulation are given in appendix A. 



III. DISCUSSION 

Figure |I] shows the invariant mass (M) distribution of unlike-sign lepton pairs in the 
ALICE detector acceptance. With a transverse momentum threshold of 2 GeV/c on each 
lepton, the dileptons from B meson decays are the dominant component of the combinatorial 
background all over the invariant mass region except for 3 ^ M ^ 5 GeV/c^ where the 
component from D meson decays is significant. Particularly interesting is the difference in 
shape between the two invariant mass distributions: the distribution from D meson decays 
exhibits one single maximum at intermediate M and a shoulder at low M while that from 
B meson decays shows two bumps at low and intermediate M. The bump at intermediate 
M results mostly from a combination of two primary leptons while the bump at low M is a 
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result of a combination of a primary and a secondary lepton produced by the same B meson 
in the so-called B-chain channel[[ This channel represents a sizeable fraction of dileptons 
because the first generation of decay products of a B meson generally contains a D meson 
whose semi-leptonic decay branching ratio is rather large. On the other hand, this channel 
is almost inexistent for primary D mesons because the secondary leptons from charm decay 
mostly originate from vr and K meson decay. Most of these secondary leptons which have 
very low momenta are removed by the transverse momentum threshold. Consequently, the 
dilepton distribution from D meson decay is dominated, to a large extent, by pairs of primary 
leptons. 

The background to the unlike-sign correlated dilepton signal is estimated by means of the 
two previously described techniques. The like-sign spectrum is normalized with the R factor 
equal to 1 since no charged particle asymmetry is present in our simulated events (more 
details about the R factor are given in appendix B) Concerning the event-mixing spectrum 
we use the standard normalization factor of 2 mentioned previously since all simulated events 
can be considered as strictly identical in terms of centrality. The resulting background 
distributions are shown by the dashed and the dotted histograms in Fig. Clearly the two 
techniques do not lead to the same result: for large M, the background from the like-sign 
method systematically underestimates that from the event-mixing method. This can be 
better observed by inspecting Fig. |^ and ^ where the background-subtracted spectrum is 
compared to the sum of all unlike-sign correlated signals. Figure § shows that the event- 
mixing technique works as expected i.e. the background-subtracted spectrum is identical 
to the sum of the unlike-sign correlated signals (the difference between the solid and the 
dashed histograms is hardly visible in Fig. |^). On the other hand. Fig. ^ shows that the 

primary lepton li is a lepton directly produced by the meson Q in the first decay generation: 
Q ^ li + anything. A secondary lepton I2 is a lepton produced in the second decay generation: 
Q — > X + anything, X ^ I2 + anything. 



6 



background-subtracted spectrum with the hke-sign technique seems to describe well the 
unlike-sign correlated signal at low M but underestimates this signal by about 25% for 
M > 3 GeV/c. 

Two effects contribute to the failure of the like-sign technique. Both of them are related 
to the particularity of the B mesons in a sense that a BB pair produces not only unlike-sign 
correlated lepton pairs, but also like-sign correlated lepton pairs: 

• The first decay generation of B mesons contains ~ 10% of primary leptons and a large 
fraction of D mesons which decay semi-leptonically with a branching ratio of ~ 12%. 
Therefore a BB pair is a source of like-sign correlated pairs. For example, in the 
following decay chain: 

B+ D° e+ z/g, — > e^ anything 
B" ^ D° TT", D° e+ anything 
the B+B^ pair produces a correlated e+e"*" pair in addition to the e+e~ pair. 

• Like the K°K° system, the 2 neutral B°B° meson systems BjB° and B^Bg undergo the 
phenomenon of particle-antiparticle mixing (or oscillation). This effect is quantified 
by the so-called mixing parameter xa (Xs) which corresponds to the time-integrated 
probability that a produced B° (B^) decays as a B° (6°) and vice versa. The mixing 
parameters are predicted by the Standard Model and have been measured experi- 
mentally |T^. They are estimated to Xd = 0.17 and Xs > 0.49. Therefore, a B^B^ 



(BgBg) pair produces, in the primary dilepton channel, ~ 70% (~ 50%) of unlike-sign 
correlated lepton pairs and ~ 30% (~ 50%) of like-sign correlated lepton pairs. 

Due to these effects the unlike-sign dilepton spectrum does not contain the full correlated 
signal because a part of this correlated signal is made of like-sign lepton pairs. Conse- 
quently, the like-sign subtraction removes from the unlike-sign dilepton spectrum not only 
the uncorrelated component but also a fraction of the correlated signal. Note that, in the 
case presented in Fig. ^, the use of the like-sign technique to subtract the background to 
the T signal would bias only weakly the result because the yield in the resonances peak is 



much larger than the yield in the background below these peaks. On the contrary, this bias 
would become important in the case of a strong T suppression which would manifest itself 
by an almost vanishing peak. It is obvious that it would become even more important for 
continuum physics. 

Figure ^ shows the like-sign and the unlike-sign correlated components. For 
M < 3 GeV/c , the correlated B meson decay distribution almost exclusively consists 
of unlike-sign pairs. For M > 3 GeV/c^ the like -sign correlated component reaches about 
45% of the unlike-sign correlated component. The like-sign correlated component arising 
from B° mixing amounts to about 30% of the total like-sign correlated distribution. 

Since D mesons do not oscillate (no evidence for D meson oscillation has been observed 
so far [|13[)) a DD pair cannot produce a like-sign correlated lepton pair in the primary decay 
generation. Nevertheless, D mesons can generate like-sign correlated lepton pairs from their 
decay chain. It can be seen from Fig. ^ that the corresponding size of the effect is very 
weak in both the dielectron channel and the dimuon channel. As already mentioned, this is 
due to the fact that most of the secondary leptons from charm decay are removed by the 
transverse momentum threshold of 2 GeV/c. 

As shown above, the event-mixing background does not contain any correlated signal 
whereas the like-sign background contains, in addition to the uncorrelated signal, the like- 
sign correlated lepton pairs from bottom decay. Therefore, subtracting the first background 
from the second one should get access to the like-sign correlated component. Figure ^ shows 
that this is indeed the case in the present simulations: when subtracting the event-mixing 
spectrum from the like-sign spectrum, one obtains a non-zero distribution which corresponds 
precisely to the like-sign correlated component from bottom decay. The latter is directly 
connected to the full bb cross-section since any change in the bb cross-section will translate 
into a proportional change in the like-sign correlated dilepton yield. The signal obtained 
from the subtraction of the two background distributions can therefore be considered as 
a reliable measurement of the number of bb pairs. Indeed, neither a resonance nor the 
thermal radiation nor the Drell-Yan mechanism can provide like-sign correlated dileptons 
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and, as shown in Fig. ^, the amount of hke-sign correlated lepton pairs from charm decay is 
neghgible when applying a pt threshold of 2 GeV/c on both leptons. 

Note that one could also get a model-independent estimate of the number of bb pairs 
from the integration of the background-subtracted unlike-sign spectrum at high invariant 
mass where the dileptons from bottom decay dominate (Fig. |I]). However, this could lead 
to a non-pure bottom sample since the integrated signal would contain all kind of unlike- 
sign correlated pairs. In particular, it would contain the Drell-Yan signal whose yield could 
become significantly large relative to the sum of correlated charm and bottom yields in the 
case of strong energy loss effects on heavy quarks . It would also contain the unlike-sign 
correlated pairs from charm whose thermal production cross-section could be enhanced for 
a QGP with a relatively high temperature |[T5|JT6|| , as well as lepton pairs from a thermalized 
hot QGP. 

We point out that whereas in e"'"e~ reactions a good understanding of open beauty hadron 
production, sample composition and decay has been recently achieved |T^, in nucleus- 
nucleus reactions none of the open beauty hadrons has ever been measured so far (this 
holds even for open charm hadrons). For the present investigations we have assumed that B 
meson production, mixing, and decay can be extrapolated from nucleon-nucleon collisions 
to nucleus-nucleus collisions by multiplication with the number of binary collisions. The 
reality may be far from this, and it is indeed a major goal of the ultra-relativistic heavy ion 
program to study the differences. However, for the present studies which concern mostly 
the signal/background we believe that the current investigation addresses the crucial and 
important new points which can be studied in the heavy ion environment. 

Note, finally, that the proposed method provides a global estimate of the like-sign signal 
from bottom including the signal from B° mixing and from the B meson decay chain. Fur- 
ther detailed investigations should reveal whether the two contributions might be identified 
separately by means of further constraints on the kinematical characteristics of the lepton 
pairs. 
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IV. CONCLUSIONS 



We have presented some features of high mass dilepton spectra in heavy ion colhsions at 
colhder energies with special emphasis on the determination of the continuum background at 
high invariant masses. Our investigations show that the combinatorial background contains 
a large amount of dileptons from bottom decay. Because bottom decay is a source of like- 
sign correlated lepton pairs, the result of the background determination using the like-sign 
technique is inappropriate. The event-mixing method is not affected by this effect and 
gives a reliable estimate of the combinatorial background. We have demonstrated that, 
by subtracting the event-mixing distribution from the like-sign distribution, one obtains a 
precise estimate of the yield of like-sign correlated lepton pairs from bottom decay. These 
results are relevant for analysis of data from RHIC and the LHC 
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APPENDIX A: SIMULATION ENVIRONMENT 



Our simulation is similar to that presented in WTj. It consists of the following steps 



Open charm and open bottom in central Pb -|- Pb reactions are computed with 



PYTHIA 5.7 [0 assuming that a Pb + Pb collision is a superposition of a certain 
number of p-|-p reactions. We use the GRVHO parton distribution function to nor- 
malize the production cross-sections in p-|-p reactions at ^/s = 5.5 TeV to 6.7 mb 
and 0.2 mb for charm and bottom respectively ||19|. These cross-sections are then 



extrapolated to Pb -|- Pb collisions at impact parameter 6 = fm with the nuclear 
overlap function TpbPb(O) = 30.4 mb^^ |2^. It results into 205 cc and 6 bb per 
central Pb + Pb event. 
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The heavy quarks hadronize into heavy mesons through a Peterson fragmentation 
function ||21|. The c(c) quarks are assumed to fragment into D+(D^), D''(D°), and 
D+(Dj7) mesons. For the b(b) quarks, we consider B+(B^), B*^(B'^), Bg(Bg), and 



Heavy mesons decay according to JETSET 6.4 with free decay branching ratios. 



Resonance events including (j), J/'4^i T, T' and T" are generated separately with 
yields and spectra taken from p2| , p3| . The resonances decay according to JETSET 
with fixed branching ratio. 

The ALICE detector response is modeled in a simple way byQ: 

— geometrical conditions: 45° < 60 < 135° and 2° < 6^ < 9° where 9e and 
6^ denote the polar angle of electrons and muons, respectively; 

— lower and upper transverse momentum thresholds: 2 < pt < 10 GeV/c for 
both electrons and muons; 

— conditions on vertex: < 3 cm and < 100 cm where vl and are the 
electron transverse distance and muon longitudinal distance between the interac- 
tion point and the track vertex. The second vertex cut simulates the effect of the 
front absorber of the muon spectrometer; 

— transverse momentum resolution: Apt/pt increases linearly from 0.6% at 
Pt=l GeV/c to 1.4% at pt=10 GeV/c for both electrons and muons. 

The leptonic component of central Pb + Pb events is constructed by means of a cocktail 
of the generated leptons from open charm, open bottom and resonances. Each lepton is 
given a weight which includes the production cross-section of its source folded with the 



^We do not intend to reproduce detailed features of the ALICE detector response but use repre- 
sentative numbers for illustrative purposes. 
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actual number of simulated sources per event and the corresponding decay branching 
ratio. The weight and the momentum components of the leptons are registered keeping 
the information on the direct parent and the grandparent of the decay product. This 
allows to trace back the composition of final spectra. 

Leptons from other sources are neglected. Not taken into account in the simulation are: 
shadowing of the structure functions, energy loss of the partons in the medium, rescattering 
of heavy mesons, suppression/enhancement of the resonance yields. The ALICE detector 
filter assumes perfect particle identification and efficiency. 

APPENDIX B: ORIGINS OF CHARGE ASYMMETRY AND THE R FACTOR 

The presence of an asymmetry in the production of positive and negative charged back- 
ground leptons leads to an enhancement of the like-sign dileptons relative to the unlike-sign 
dileptons. Therefore, the correct normalization of the like-sign spectrum relies on a precise 
determination of the R factor. We discuss here the origins of possible charge asymmetries in 
nucleus-nucleus reactions, in which conditions this leads to a R factor different from unity 
and how to estimate R in this case. 

Charge asymmetry at the production level 

At least three physical effects at the particle production level can lead to a charge asym- 
metry in the final stage of nucleus-nucleus reactions: 

• charge conservation implies an excess of positively charged hadrons in the final state 
because the colliding ions are positively charged; 

• isospin conservation implies a ratio tt^/tt'*' > 1 because the ion N/Z ratio is larger 
than 1; 

• associate Kaon production results in an excess of K+ relative to K~. 
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It is obvious that these asymmetries propagate to the decayed leptons. However in central 
nucleus-nucleus reactions the multiplicity of charged particles produced in the interaction 
is so large that any initial charge asymmetry is smeared-out In fact, it can be demon- 
strated that R = 1 exactly if the particle multiplicities are Poisonnian This has also 
been shown by means of Monte-Carlo simulations for central Pb + Pb reactions at SPS |^ 



The charge asymmetry tends to vanish with increasing beam energy as evidenced by the mea- 
sured tt^/tt^ ratio which goes, in central Au + Au reactions, from 2.15±0.30 at GSI/SIS ||2^ 
to 1.00±0.02 at RHIC 

On the contrary, the previous statement is not valid anymore in low multiplicity events 
like non-central nucleus-nucleus reactions. In this case R is larger than unity. Its actual value 
can be precisely determined by means of simulations [^. Even in this case R differs only 
slightly from unity. Indeed in peripheral (b ~ 13 fm) Pb + Pb collisions at SPS, the R factor 
for dimuons in the acceptance of the NA50 spectrometer is found to be 1.075 |2^. It should 



be even closer to unity at higher beam energy thanks to the larger particle multiplicity. 

Charge asymmetry due to a detector acceptance/trigger bias 

A charge asymmetry in the lepton sample can also be the consequence of a possible 
detector bias relative to the lepton charge. The detector bias caused by a different acceptance 
for positively and negatively charged particle is usually taken into account by averaging the 
data collected with the two opposite settings of the magnetic field p5|. Similarly a charge 
asymmetry could result from a different probability for the trigger system to accept like-sign 
and unlike-sign events. Such an effect can be estimated by means of simulations as discussed 
in Ml. 
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FIG. 1. Invariant mass distributions of unlike-sign dielectrons (left) and unlike-sign dimuons 
(right) for central Pb + Pb reactions at LHC energy. The events are filtered through the ALICE de- 
tector acceptance cuts as described in the text. The top (solid) histogram on each figure comprises 
the sum of all unlike-sign lepton pairs. The dashed and dotted histograms show the background 
determined by the event-mixing technique and with the like-sign technique, respectively. The his- 
tograms labeled charm and bottom correspond to the component where both leptons result from 
D and B meson decay, respectively. The two lower panels are zooms of the upper panels in the 
invariant mass region of the T family. 
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FIG. 2. Invariant mass distributions of unlike-sign correlated dileptons (solid histograms) for 
central Pb + Pb reactions at LHC energy. The events are filtered through the ALICE detector 
acceptance cuts as described in the text. The dashed histogram is obtained by the subtraction of 
the background, determined by the event-mixing technique, from the total distribution. For more 
details, see text. 
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FIG. 3. Same as Fig. ^ but using the like-sign technique instead of the event-mixing technique 
for the determination of the background. 
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FIG. 4. Invariant mass distributions of correlated dileptons from bottom decay for central 
Pb + Pb reactions at LHC energy. The events are filtered through the ALICE detector acceptance 
cuts as described in the text. The unlike-sign and like-sign components are shown by the solid and 
dashed histograms, respectively. The dotted histogram represents the like-sign component which 
results from B'^B'' mixing. 
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FIG. 5. Invariant mass distributions of correlated dileptons from charm decay for central 
Pb + Pb reactions at LHC energy. The events are filtered through the ALICE detector accep- 
tance cuts as described in the text. The unlike-sign and like-sign components are shown by the 
solid and dashed histograms, respectively. 
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FIG. 6. Invariant mass distributions of dileptons for central Pb + Pb reactions at LHC energy. 
The events are filtered through the ALICE detector acceptance cuts as described in the text. Solid 
histogram: like-sign correlated dileptons from B meson decays. Dashed histogram: distributions 
obtained after subtraction of the background estimated with the event-mixing technique from the 
one estimated by the like-sign technique. 
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